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A consistent physical theory of quantum mechanics can be built on a complex 
Hamiltonian that is not Hermitian but instead satisfies the physical condition of 
space-time reflection symmetry (VT symmetry). Thus, there are infinitely many 
new Hamiltonians that one can construct that might explain experimental data. One 
would think that a quantum theory based on a non-Hermitian Hamiltonian violates 
unitarity. However, if VT symmetry is not broken, it is possible to use a previously 
unnoticed physical symmetry of the Hamiltonian to construct an inner product whose 
associated norm is positive definite. This construction is general and works for any 
"PT-symmetric Hamiltonian. The dynamics is governed by unitary time evolution. 
This formulation does not conflict with the requirements of conventional quantum 
mechanics. There are many possible observable and experimental consequences of 
extending quantum mechanics into the complex domain, both in particle physics and 
in solid state physics. 
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I. INTRODUCTION 

In this paper we present an alternative to the one of the standard axioms of quantum 
mechanics; namely, that the Hamiltonian H, which incorporates the symmetries and specifies 
the dynamics of a quantum theory, must be Hermitian: H = W. It is commonly believed 
that the Hamiltonian must be Hermitian in order to ensure that the energy spectrum (the 
eigenvalues of the Hamiltonian) is real and that the time evolution of the theory is unitary 
(probability is conserved in time). Although this axiom is sufficient to guarantee these 
desired properties, we argue here that it is not necessary. We believe that the condition 
of Hermiticity is a mathematical requirement whose physical basis is somewhat remote and 
obscure. We demonstrate here that there is a simpler and more physical alternative axiom, 
which we refer to as space-time reflection symmetry (VT symmetry): H = H vr . This 
symmetry allows for the possibility of non-Hermitian and complex Hamiltonians but still 
leads to a consistent theory of quantum mechanics. 

We also show that because VT symmetry is an alternative condition to Hermiticity 
it is now possible to construct infinitely many new Hamiltonians that would have been 
rejected in the past because they are not Hermitian. An example of such a Hamiltonian is 
H = p 2 + ix 3 . It should be emphasized that we do not regard the condition of Hermiticity 
as wrong. Rather, the condition of VT symmetry offers the possibility of studying new and 
interesting quantum theories. 

Let us recall the properties of the space reflection (parity) operator V and the time- 
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reflection operator T . The parity operator V is linear and has the effect 

p — > — p and x — * —x. 

The time-reversal operator T is antilinear and has the effect 

p — > —p, x — > x, and i — > —i. 

Note that T changes the sign of i because, like the parity operator, it preserves the fun- 
damental commutation relation of quantum mechanics, [x,p] = i, known as the Heisenberg 
algebra |26j]. 

It is easy to construct infinitely many Hamiltonians that are not Hermitian but do possess 
VT symmetry. For example, consider the one-parameter family of Hamiltonians 

H = p 2 + x 2 {ix) e (ereal). (1) 

Note that while H in is not symmetric under V or T separately, it is invariant under their 
combined operation. We say that such Hamiltonians possess space-time reflection symmetry. 
Other examples of complex Hamiltonians having VT symmetry are H = p 2 + x A (ix) e , 
H = p 2 + x 6 (ix) e , and so on 00. 

The class of PT-symmetric Hamiltonians is larger than and includes real symmetric 
Hermitians because any real symmetric Hamiltonian is automatically PT-symmetric. For 
example, consider the real symmetric Hamiltonian H = p 2 + x 2 + 2x. This Hamiltonian is 
time-reversal symmetric, but according to the usual definition of space reflection for which 
x — ► —x, this Hamiltonian appears not to have VT symmetry. However, recall that the 
parity operator is defined only up to unitary equivalence 0. In this example, if we express 
the Hamiltonian in the form H = p 2 + (x + l) 2 — 1, then it is evident that H is VT symmetric, 
provided that the parity operator performs a space reflection about the point x — — 1 rather 
than x = 0. See Ref. [l| for the general construction of the relevant parity operator. 

Five years ago it was discovered that with properly defined boundary conditions the 
spectrum of the Hamiltonian H in (JTJ is real and positive when e > The spectrum is 
partly real and partly complex when e < 0. The eigenvalues have been computed numerically 
to very high precision, and the real eigenvalues are plotted as functions of e in Fig. Q 

We say that the VT symmetry of a Hamiltonian H is unbroken if all of the eigenfunctions 
of H are simultaneously eigenfunctions of VT [28]. It is easy to show that if the VT 
symmetry of a Hamiltonian H is unbroken, then the spectrum of H is real. The proof is 
short and goes as follows: Assume that a Hamiltonian H possesses VT symmetry (that is, 
that H commutes with the VT operator), and that if <fi is an eigenstate of H with eigenvalue 
E, then it is simultaneously an eigenstate of VT with eigenvalue A: 

H<j> = E(f) and VT = \<f>. (2) 

We begin by showing that the eigenvalue A is a pure phase. Multiplying VTip = A0 on 
the left by VT and using the fact that V and T commute and that V 2 = T 2 = 1 we conclude 
that = A*A0 and thus A = e lOL for some real a. Next, we introduce the convention that 
is used throughout this paper. Without loss of generality we replace the eigenstate by 
e -ia/20 go faofc it s eigenvalue under the operator VT is unity: 



VT<fi = <f). 



(3) 
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FIG. 1: Energy levels of the Hamiltonian H = p 2 + x 2 (ix) e as a function of the parameter e. There 
are three regions: When e > 0, the spectrum is real and positive and the energy levels rise with 
increasing e. The lower bound of this region, e = 0, corresponds to the harmonic oscillator, whose 
energy levels are E n = 2n + 1. When — 1 < e < 0, there are a finite number of real positive 
eigenvalues and an infinite number of complex conjugate pairs of eigenvalues. As e decreases from 
to —1, the number of real eigenvalues decreases; when e < —0.57793, the only real eigenvalue 
is the ground-state energy. As e approaches — 1 + , the ground-state energy diverges. For e < — 1 
there are no real eigenvalues. 

Let us turn to the eigenvalue equation Hep = E<p. We multiply this equation on the left 
by VT and use the fact that [VT, H] = to obtain E<p = E*<p. Hence, E = E* and the 
eigenvalue E is real. 

The crucial assumption in this argument is that is simultaneously an eigenstate of H 
and VT. In quantum mechanics if a linear operator X commutes with the Hamiltonian 
H, then the eigenstates of H are also eigenstates of X. However, we emphasize that the 
operator VT is not linear (it is antilinear) and thus we must make the extra assumption 
that the VT symmetry of H is unbroken; that is, that <fi is simultaneously an eigenstate of 
H and VT. This extra assumption is nontrivial because it is not easy to determine a priori 
whether the VT symmetry of a particular Hamiltonian H is broken or unbroken. For the 
Hamiltonian H in (0) the VT symmetry is unbroken when e > and it is broken when e < 0. 
Note that the conventional Hermitian Hamiltonian for the quantum mechanical harmonic 
oscillator lies at the boundary of the unbroken and the broken regimes. Recently, Dorey et 
al. proved rigorously that the spectrum of H in is real and positive (Hf in the region 
e > 0. Many other PT-symmetric Hamiltonians for which space-time reflection symmetry 
is not broken have been investigated, and the spectra of these Hamiltonians have also been 
shown to be real and positive [6|. 

While it is useful to show that a given non- Hermitian PT-symmetric Hamiltonian opera- 
tor has a positive real spectrum, the urgent question that must be answered is whether such 
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a Hamiltonian defines a physical theory of quantum mechanics. By a physical theory we 
mean that there is a Hilbert space of state vectors and that this Hilbert space has an inner 
product with a positive norm. In the theory of quantum mechanics we interpret the norm of 
a state as a probability and this probability must be positive. Furthermore, we must show 
that the time evolution of the theory is unitary. This means that as a state vector evolves 
in time the probability does not leak away. 

It is not at all obvious whether a Hamiltonian such as H in (P) gives rise to a consistent 
quantum theory. Indeed, while past investigations of this Hamiltonian have shown that the 
spectrum is entirely real and positive when e > 0, it appeared that one inevitably encoun- 
tered the severe problem of dealing with Hilbert spaces endowed with indefinite metrics 0. 
In this paper we will identify a new symmetry that all PT-symmetric Hamiltonians having 
an unbroken PT-symmetry possess. We denote the operator representing this symmetry by 
C because the properties of this operator resemble those of the charge conjugation operator 
in particle physics. This will allow us to introduce an inner product structure associated 
with CVT conjugation for which the norms of quantum states are positive definite. We 
will see that CVT symmetry is an alternative to the conventional Hermiticity requirement; 
it introduces the new concept of a dynamically determined inner product (one that is de- 
fined by the Hamiltonian itself). As a consequence, we will extend the Hamiltonian and 
its eigenstates into the complex domain so that the associated eigenvalues are real and the 
underlying dynamics is unitary. 



II. CONSTRUCTION OF THE C OPERATOR 

We begin by summarizing the mathematical properties of the solution to the Sturm- 
Liouville differential equation eigenvalue problem 

-<j)'n{x) + x 2 (ix) e 4> n (x) = E n (p n (x) (4) 

associated with the Hamiltonian H in (0). The differential equation (0J) must be imposed 
on an infinite contour in the complex-x plane. For large \x\ this contour lies in wedges that 
are placed symmetrically with respect to the imaginary-x axis 3] . The boundary conditions 
on the eigenf unctions are that <j)(x) — *> exponentially rapidly as |x| — > oo on the contour. 
For < e < 2, the contour may be taken to be the real axis. 

When e > 0, the Hamiltonian has an unbroken VT symmetry. Thus, the eigenf unctions 
4> n (x) are simultaneously eigenstates of the VT operator: VT<p n (x) = \ n (p n (x). As we 
argued above, A n is a pure phase and, without loss of generality, for each n this phase can 
be absorbed into (p n (x) by a multiplicative rescaling so that the new eigenvalue is unity: 

VT<\> n {x) = <f>* n {-x) = <f> n {x). (5) 

There is strong evidence that, when properly normalized, the eigenfunctions <p n {x) are 
complete. The coordinate-space statement of completeness (for real x and y) reads 

J2(-i) n M^)My) = 5(x- y ). (6) 

This is a nontrivial result that has been verified numerically to extremely high accuracy 
(twenty decimal places) jHEJ. Note that there is a factor of (— l) n in the sum. This unusual 
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factor does not appear in conventional quantum mechanics. The presence of this factor is 
explained in the following discussion of ortho normality [see (JHJ)]. 

Here is where we encounter the underlying problem associated with non-Hermitian VT- 
symmetric Hamiltonians. There seems to be a natural choice for the inner product of two 
functions f(x) and g(x): 

(f,g) = J dx[PTf(x)]g(x), (7) 

where VTf(x) = [/(— x)]* and the integral is taken over the contour described above in 
the complex-x plane. The apparent advantage of this inner product is that the associated 
norm (/, /) is independent of the overall phase of f(x) and is conserved in time. Phase 
independence is desired because in the theory of quantum mechanics the objective is to 
construct a space of rays to represent quantum mechanical states. With respect to this 
inner product the eigenfunctions (j) m {x) and <p n {x) of H in {1} are orthogonal for n ^ m. 
However, when m = n the norm is evidently not positive: 

(<PmAn) = (-l) n Smn. (8) 

This result is apparently true for all values of e in (J3J) and it has been verified numerically 
to extremely high precision. Because the norms of the eigenfunctions alternate in sign, 
the Hilbert space metric associated with the VT inner product (•, •) is indefinite. This 
split signature (sign alternation) is a generic feature of the VT inner product. Extensive 
numerical calculations verify that the formula in (jHJ) holds for all e > 0. 

Despite the lack of positivity of the inner product, we proceed with the usual analysis 
that one would perform for any Sturm-Liouville problem of the form H(f) n = E n <p n . First, 
we use the inner product formula (JHJ) to verify that (JBJ) is the representation of the unity 
operator. That is, we verify that 

dyS(x — y)6(y — z) — S(x — z). (9) 

Second, we reconstruct the parity operator V in terms of the eigenstates. The parity operator 
in position space is V(x, y) = 5(x + y), so from (JHJ) we get 

v(x, y) = J2(-i) n Mx)M-y)- (io) 

n 

By virtue of (JHJ) the square of the parity operator is unity: V 2 = 1. 
Third, we reconstruct the Hamiltonian H in coordinate space: 

H(x,y) = Y,(-l) n E n Mx)My)- (11) 

n 

Using (JUJ) - (JHJ) it is easy to see that this Hamiltonian satisfies H(j) n (x) = E n <p n (x). Fourth, 
we construct the coordinate-space Green's function G(x,y): 

G(x, y) = J2(-l) n -^rMx)My)- (12) 

Note that the Green's function is the functional inverse of the Hamiltonian; that is, G 
satisfies the equation 

d 2 



dyH(x,y)G(y,z) 



dx 2 



+ x 2 (ix) 



G{x,z) = 5{x - z). (13) 
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While the time-independent Schrodinger equation (jlj) cannot be solved analytically, the 
differential equation for G(x, z) in (JT3j) can be solved exactly and in closed form 0. The 
technique is to consider the case < e < 2 so that we may treat x as real and then to 
decompose the x axis into two regions, x > z and x < z. We can solve the differential 
equation in each of these regions in terms of Bessel functions. Then, using this coordinate- 
space representation of the Green's function, we construct an exact closed-form expression 
for the spectral zeta function (sum of the inverses of the energy eigenvalues). To do so we 
set y — x in G(x, y) and use (jSJ) to integrate over x. For all e > we obtain 



E- 

n 



cos (jfy sin (£) 



(* + ')»r(g)r(g) 



Having presented these general Sturm-Liouville constructions, we now address the crucial 
question of whether a PT-symmetric Hamiltonian defines a physically viable quantum me- 
chanics or whether it merely provides an intriguing Sturm-Liouville eigenvalue problem. The 
apparent difficulty with formulating a quantum theory is that the vector space of quantum 
states is spanned by energy eigenstates, of which half have norm +1 and half have norm —1. 
Because the norm of the states carries a probabilistic interpretation in standard quantum 
theory, the existence of an indefinite metric in (jSJ) seems to be a serious obstacle. 

The situation here in which half of the energy eigenstates have positive norm and half 
have negative norm is analogous to the problem that Dirac encountered in formulating the 
spinor wave equation in relativistic quantum theory |l0|. Following Dirac's approach, we 
attack the problem of an indefinite norm by finding a physical interpretation for the negative 
norm states. We claim that in any theory having an unbroken VT symmetry there exists 
a symmetry of the Hamiltonian connected with the fact that there are equal numbers of 
positive-norm and negative-norm states. To describe this symmetry we construct a linear 
operator denoted by C and represented in position space as a sum over the energy eigenstates 
of the Hamiltonian [111 ]: 

C(x,y) = ^(j) n (x)(j) n (y). (15) 

n 

As stated earlier, the properties of this new operator C are nearly identical to those of 
the charge conjugation operator in quantum field theory. For example, we can use equations 
© - © to verify that the square of C is unity (C 2 = 1): 

J dyC(x,y)C(y,z) = 5(x - z). (16) 

Thus, the eigenvalues of C are ±1. Also, C commutes with the Hamiltonian H. Therefore, 
since C is linear, the eigenstates of H have definite values of C. Specifically, if the energy 
eigenstates satisfy (jHJ), then we have C(j) n = (—l) n (f> n because 



C(p n (x) = J dyC(x,y)(j) n {y) = J^. 



4>m(x) / dy<f) m (y)<f) n (y). 



We then use J dy (p m (y)<p n (y) = (<f) m , <fi n ) according to our convention. We conclude that C 
is the operator observable that represents the measurement of the signature of the VT norm 
of a state l2£ 
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Note that the operators V and C are distinct square roots of the unity operator 8{x — y). 
That is, while V 2 = 1 and C 2 = 1, V and C are not identical. Indeed, the parity operator V 
is real, while C is complex [30]. Furthermore, these two operators do not commute; in the 
position representation 

{CP)(x,y) = J2<f>n(x)<f>n{-y) but (VC)(x,y) = J2M-x)My), (17) 

n n 

which shows that CV = (VC)*. However, C does commute with VT. 

Finally, having obtained the operator C we define a new inner product structure having 
positive definite signature by 

(f\g}= [ dx[CVTf(x)}g(x). (18) 

Like the VT inner product (J7J), this inner product is phase independent and conserved in 
time. This is because the time evolution operator, just as in ordinary quantum mechanics, 
is e lHt . The fact that H commutes with the VT and the CVT operators implies that both 
inner products, (JJJ) and (|18|). remain time independent as the states evolve in time. However, 
unlike (JZJ), the inner product (fT8j) is positive definite because C contributes —1 when it acts on 
states with negative VT norm. In terms of the CVT conjugate, the completeness condition 
flU) reads 

Y,<Pn{x)[CVT<p n {y)} =6(x-y). (19) 

n 

Unlike the inner product of conventional quantum mechanics, the CVT inner product (|19|) 
is dynamically determined; it depends implicitly on the choice of Hamiltonian. 

The operator C does not exist as a distinct entity in conventional quantum mechanics. 
Indeed, if we allow the parameter e in (|TJ) to tend to zero, the operator C in this limit 
becomes identical to V. Thus, in this limit the CVT operator becomes T, which is just 
complex conjugation. As a consequence, the inner product (fTHj) defined with respect to the 
CVT conjugation reduces to the complex conjugate inner product of conventional quantum 
mechanics when e — > 0. Similarly, in this limit (|T9~|) reduces to the usual statement of 
completeness Y^ n ^{x)(jf n {y) = 6(x - y). 

Note that the CVT inner-product (fTHjl is independent of the choice of integration contour 
C so long as C lies inside the asymptotic wedges associated with the boundary conditions for 
the Sturm-Liouville problem (J2J). Path independence is a consequence of Cauchy's theorem 
and the analyticity of the integrand. In conventional quantum mechanics, where the positive- 
definite inner product has the form J dx f*(x)g(x), the integral must be taken along the real 
axis and the path of the i nteg ration cannot be deformed into the complex plane because the 



integrand is not analytic j3l|. The VT inner product (J7J) shares with (fTHj) the advantage of 
analyticity and path independence, but suffers from nonpositivity. We find it surprising that 
a positive-definite metric can be constructed using CVT conjugation without disturbing the 
path independence of the inner-product integral. 

Finally, we explain why PT-symmetric theories are unitary. Time evolution is deter- 
mined by the operator e~ lHt , whether the theory is expressed in terms of a PT-symmetric 
Hamiltonian or just an ordinary Hermitian Hamiltonian. To establish the global unitarity 
of a theory we must show that as a state vector evolves its norm does not change in time. 
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If ipo(x) is a prescribed initial wave function belonging to the Hilbert space spanned by the 
energy eigenstates, then it evolves into the state ipt( x ) at time t according to 

VtO) = e- iHt ip (x). 

With respect to the CVT inner product defined in (j!8|) . the norm of the vector ipt{ x ) does 
not change in time, 

hl>t\*l>t) = (V'ol^o), 

because the Hamiltonian H commutes with the CVT operator. Establishing unitarity at 
a local level is more difficult. Here, we must show that in coordinate space, there exists a 
local probability density that satisfies a continuity equation so that the probability does not 
leak away. This is a subtle result because the probability current flows about in the complex 
plane rather than along the real axis as in conventional Hermitian quantum mechanics. 
Preliminary numerical studies indeed indicate that the continuity equation is fulfilled [liij . 



III. ILLUSTRATIVE EXAMPLE: A 2 x 2 MATRIX HAMILTONIAN 

We will now illustrate the above results concerning PT-symmetric quantum mechanics in 
a very simple context. To do so we will consider systems characterized by finite-dimensional 
matrix Hamiltonians. In finite-dimensional systems the V, T, and C operators appear, but 
there is no analogue of the boundary conditions associated with coordinate-space Schrodinger 
equations. 

Let us consider the 2x2 matrix Hamiltonian 

*=( r f <»> 

where the three parameters r, s, and 9 are real. This Hamiltonian is not Hermitian in the 
usual sense, but it is VT symmetric, where the parity operator is given by [lij 

and T performs complex conjugation. 

There are two parametric regions for this Hamiltonian. When s 2 < r 2 sm 2 9, the energy 
eigenvalues form a complex conjugate pair. This is the region of broken VT symmetry. On 
the other hand, if s 2 > r 2 sin 2 9, then the eigenvalues e± = r cos 9±^ — r 2 sin 9 are real. 
This is the region of unbroken VT symmetry. In the unbroken region the simultaneous 
eigenstates of the operators H and VT are given by 

*+> = f and |e_) = ~ 7 == * , (22) 



V2^\e~ ia/ V 1 ' v/2 



cos a 



where we set sin a = (r/s) sin 9. It is easily verified that (e±, e±) = ±1 and that (e±, e T ) = 0, 
recalling that (u, v) = {VTu) ■ v. Therefore, with respect to the VT inner product, the 
resulting vector space spanned by energy eigenstates has a metric of signature (+, — ). The 
condition s 2 > r 2 sin 2 9 ensures that VT symmetry is not broken. If this condition is violated, 
the states (|2*2^1 are no longer eigenstates of VT because a becomes imaginary [32|. 
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Next, we construct the operator C: 

C= — ( i ^ a 1 Y (23) 
cos a V 1 sin ay 

Note that C is distinct from H and V and has the key property that 

C|e±> = ±|e±>. (24) 

The operator C commutes with H and satisfies C 2 = 1. The eigenvalues of C are precisely 
the signs of the VT norms of the corresponding eigenstates. 

Using the operator C we construct the new inner product structure 

(u\v) = (CVTu) -v. (25) 

This inner product is positive definite because (e±|e±) = 1. Thus, the two-dimensional 
Hilbert space spanned by \e±), with inner product (-|-), has a Hermitian structure with 
signature (+, +). 

Let us demonstrate explicitly that the CVT norm of any vector is positive. We choose 
the arbitrary vector ip = (®), where a and b are any complex numbers. We then see 
that = (£), that VT^ = (£), and that CVT^ = ^ Thus, <^> = 

(CVTtp) ■ if) = + b*b + i(b*b — a* a) sin a]. Now let a = x + iy and b = u + iv, where 

x, y, u, and v are real. Then 

(ip\ip) = (x 2 + v 2 + 2xv sin a + y 2 + u 2 — 2yu sin a) , (26) 

cos a v ' 

which is explicitly positive and vanishes only ifx = y = u = v = 0. 

Recalling that (u\ denotes the CPT-conjugate of \u), the completeness condition reads 

\e + )(e + \ + \eJ)(e-.\= jV (27) 

Furthermore, using the CVT conjugate (e±|, we can express C in the form C = \s+)(e+\ — 
|e_)(e_|, as opposed to the representation in (JT3J), which uses the VT conjugate. 

In general, an observable in this theory is represented by a CVT invariant operator; that 
is, one that commutes with CVT . Thus, if CVT symmetry is not broken, the eigenvalues 
of the observable are real. The operator C satisfies this requirement, and hence it is an 
observable. For the two-state system, if we set 6 = 0, then the Hamiltonian (f2*U|) becomes 
Hermitian. However, the operator C then reduces to the parity operator V. As a conse- 
quence, the requirement of CVT invariance reduces to the standard condition of Hermiticity 
for a symmetric matrix, namely, that H = H*. This is why the hidden symmetry C was 
not noticed previously. The operator C emerges only when we extend a real symmetric 
Hamiltonian into the complex domain. 

We have also calculated the C operator in infinite-dimensional quantum mechanical mod- 
els. For an x 2 + ix 3 potential C can be obtained from the summation in (|15J) using per- 
turbative methods and for an x 2 — x 4 potential C can be calculated using nonperturbative 
methods 
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IV. APPLICATIONS AND POSSIBLE OBSERVABLE CONSEQUENCES 

In summary, we have described an alternative to the axiom of Hermiticity in quantum 
mechanics; we call this new requirement VT invariance. In quantum field theory, Hermitic- 
ity, Lorentz invariance, and a positive spectrum are crucial for establishing CVT invariance 
Here, we have established the converse of the CVT theorem in the following limited 
sense: We assume that the Hamiltonian possesses space-time reflection symmetry, and that 
this symmetry is not broken. From these assumptions, we know that the spectrum is real 
and positive and we construct an operator C that is like the charge conjugation operator. 
We show that quantum states in this theory have positive norms with respect to CVT 
conjugation. In effect, we replace the mathematical condition of Hermiticity, whose phys- 
ical content is somewhat remote and obscure, by the physical condition of space-time and 
charge-conjugation symmetry. These symmetries ensure the reality of the spectrum of the 
Hamiltonian in complex quantum theories. 

Could non-Hermitian, PT-symmetric Hamiltonians be used to describe experimentally 
observable phenomena? Non-Hermitian Hamiltonians have already been used to describe 
interacting systems. For example in 1959, Wu showed that the ground state of a Bose 
system of hard spheres is described by a non-Hermitian Hamiltonian ^6|. Wu found that 
the ground-state energy of this system is real and conjectured that all of the energy levels 
were real. In 1992, Hollowood showed that even th oug h the Hamiltonian of a complex Toda 
lattice is non-Hermitian, the energy levels are real [17[ • Non-Hermitian Hamiltonians of the 
form H = p 2 + ix 3 also arise in various Reggeon field theory models that exhibit real positive 
spectra [18]. In each of these cases the fact that a non-Hermitian Hamiltonian had a real 
spectrum appeared mysterious at the time, but now the explanation is simple: In each of 
these cases it is easy to show that the non-Hermitian Hamiltonian is PT-symmetric. That 
is, the Hamiltonian in each case is constructed so that the position operator x or the field 
operator is always multiplied by i. 

An experimental signal of a complex Hamiltonian might be found in the context of con- 
densed matter physics. Consider the complex crystal lattice whose potential is given by 
V(x) = i sin x. While the Hamiltonian H = p 2 +i sin x is not Hermitian, it is PT-symmetric, 
and all of the energy bands are real. However, at the edge of the bands the wave function of 
a particle in such a lattice is always bosonic (27r-periodic) and, unlike the case of ordinary 
crystal lattices, the wave function is never fermionic (47r-periodic) Direct observation 
of such a band structure would give unambiguous evidence of a PT-symmetric Hamiltonian. 

There are many opportunities for the use of non-Hermitian Hamiltonians in the study 
of quantum field theory. For example, a scalar quantum field theory with a cubic self- 
interaction described by the Lagrangian C = |(V<£>) 2 + \m 2 ip 2 + gtp 3 is physically unac- 
ceptable because the energy spectrum is not bounded below. However, the cubic scalar 
quantum field theory that corresponds to H in ([Tj) with e = 1 is given by the Lagrangian 
density C = |(Vv?) 2 + \m 2 ip 2 + igip 3 . This is a new, physically acceptable quantum field 
theory. Moreover, the theory that corresponds to H in (JH) with e = 2 is described by the 
Lagrangian density 

£=^) 2 + imV-^ 4 (28) 

What is remarkable about this "wrong-sign" field theory is that, in addition to the energy 
spectrum being real and positive, the one-point Green's function (the vacuum expectation 
value of the field ip) is nonzero |2i|. Furthermore, the field theory is renormalizable, and in 
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four dimensions is asymptotically free (and thus nontrivial) [2l| . Based on these features of 
the theory, we believe that the theory may provide a useful setting to describe the dynamics 
of the Higgs sector in the standard model. 

Other field theory models whose Hamiltonians are non-Hermitian and PT-symmetric 
have also been studied. For example, PT-symmetric electrodynamics is particularly inter- 
esting because it is asymptotically free (unlike ordinary electrodynamics) and because the 
direction of the Casimir force is the negative of that in ordinary electrodynamics (13] . This 
theory is remarkable because it can determine its own coupling constant. Supersymmetric 
'PT-symmetric quantum field theories have also been studied [23(. 

We have found that PT-symmetric quantum theories exhibit surprising and new phenom- 
ena. For example, when g is sufficiently small, the —g^p A theory described by the Lagrangian 
fl2~%j) possesses bound states (the conventional g(p 4 theory does not because the potential is 
repulsive). The bound states occur for all dimensions < D < 3 j2J], but for purposes 
of illustration we describe the bound states in the context of one-dimensional quantum 
field theory (quantum mechanics). For the conventional quantum mechanical anharmonic 
oscillator, which is described by the Hamiltonian 

H = l -p 2 + X -m 2 x 2 + ^gx 4 {g > 0), (29) 

the small-g Rayleigh-Schrodinger perturbation series for the kth energy level is 

1 3 



k + ~ + ^(2k 2 + 2k + \)v + 0(z/ 2 ) 



+ ), (30) 



where v = g/(4m 3 ). The renormalized mass M is defined as the first excitation above the 
ground state: M = E x - E ~ m[l + 3v + 0(z/ 2 )] as v -> 0+ 

To determine if the two-particle state is bound, we examine the second excitation above 
the ground state using (j3~Uj). We define 

B 2 = E 2 - E ~ m [2 + 9v + 0(u 2 )] {y -> + ). (31) 

If B 2 < 2M, then a two-particle bound state exists and the (negative) binding energy is 
B 2 — 2M. If B 2 > 2M, then the second excitation above the vacuum is interpreted as 
an unbound two-particle state. We see from (|31|) that in the small-coupling region, where 
perturbation theory is valid, the conventional anharmonic oscillator does not possess a bound 
state. Indeed, using WKB, variational methods, or numerical calculations, one can show 
that there is no two-particle bound state for any value of g > 0. Because there is no bound 
state the gx A interaction may be considered to represent a repulsive force |33j|. 

We obtain the perturbation series for the non-Hermitian, PT-symmetric Hamiltonian 

H = l -p 2 + l -m 2 x 2 - -^gx 4 (g > 0), (32) 

from the perturbation series for the conventional anharmonic oscillator by replacing v with 
—v. Thus, while the conventional anharmonic oscillator does not possess a two-particle 
bound state, the PT-symmetric oscillator does indeed possess such a state. We measure 
the binding energy of this state in units of the renormalized mass M and we define the 
dimensionless binding energy A 2 by 

A 2 = B2 M 2M ~ -31/ + Q(^ 2 ) (z/^0+). (33) 
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This bound state disappears when v increases beyond v = 0.0465 .... As v continues to 
increase, A 2 reaches a maximum value of 0.427 at v = 0.13 and then approaches the limiting 
value 0.28 as v — > oo. 

In the PT-symmetric anharmonic oscillator, there are not only two-particle bound states 
for small coupling constant but also A;-particle bound states for all k > 2. The dimensionless 
binding energies are 

A k = (B k - kM)/M ~ -3k(k - l)v/2 + 0(z/ 2 ) (i/->0+). (34) 

The key feature of this equation is that the coefficient of v is negative. Since the dimension- 
less binding energy becomes negative as v increases from 0, there is a /c-particle bound state. 
The higher multiparticle bound states cease to be bound for smaller values of v\ starting 
with the three-particle bound state, the binding energy of these states becomes positive as 
v increases past 0.039, 0.034, 0.030, and 0.027. 

Thus, for any value of v there are always a finite number of bound states and an infinite 
number of unbound states. The number of bound states decreases with increasing v until 
there are no bound states at all. There is a range of v for which there are only two- and 
three-particle bound states. This situation is analogous to the physical world in which 
one observes only states of two and three bound quarks. In this range of v if one has 
an initial state containing a number of particles (renormalized masses), these particles will 
clump together into bound states, releasing energy in the process. Depending on the value 
of is, the final state will consist either of two- or of three-particle bound states, whichever 
is energetically favored. Note also that there is a special value of v for which two- and 
three-particle bound states can exist in thermodynamic equilibrium. 

How does a gy? 3 theory compare with a gip 4 theory? A gtp 3 theory has an attractive force. 
The bound states that arise as a consequence of this force can be found by using the Bethe- 
Salpeter equation. However, the gy? 3 field theory is unacceptable because the spectrum is 
not bounded below. If we replace g by ig, the spectrum becomes real and positive, but now 
the force becomes repulsive and there are no bound states. The same is true for a two-scalar 
theory with interaction of the form ig<p 2 X- This latter theory is an acceptable model of 
scalar electrodynamics, but has no analog of positronium. 

Another feature of PT-symmetric quantum field theory that distinguishes it from the 
conventional quantum field theory lies in the commutation relation between the V and C 
operators. Specifically, if we write C = Cr + iC\, where Cr and C\ are real, then CrP = VCr 
and C{P = —VC\. These commutation and anticommutation relations suggest the possi- 
bility of interpreting PT-symmetric quantum field theory as describing both bosonic and 
fermionic degrees of freedom, an idea analogous to the supersymmetric quantum theories. 
The distinction here, however, is that the supersymmetry can be broken; that is, bosonic 
and fermionic counterparts can have different masses without breaking the VT symmetry. 
Therefore, another possible observable experimental consequence might be the breaking of 
the supersymmetry. 

V. CONCLUDING REMARKS 

We have argued in this paper that there is an alternative to the axiom of standard 
quantum mechanics that the Hamiltonian must be Hermitian. We have shown that the axiom 
of Hermiticity may be replaced by the more physical condition of VT (space-time reflection) 
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symmetry. Space-time reflection symmetry is distinct from the condition of Hermiticity, so it 
is possible to consider new kinds of quantum theories, such as quantum field theories whose 
self-interaction potentials are igip 3 or —gtp A . Such theories have previously been thought to 
be mathematically and physically unacceptable because the spectrum might not be real and 
because the time evolution might not be unitary. 

These new kinds of theories can be thought of as extensions of ordinary quantum mechan- 
ics into the complex plane; that is, continuations of real symmetric Hamiltonians to complex 
Hamiltonians. The idea of analytically continuing a Hamiltonian was first discussed in 1952 
by Dyson, who argued heuristically that perturbation theory for quantum electrodynamics 
is divergent |25| . Dyson's argument involves rotating the electric charge e into the complex 
plane e — > ie. Applied to the quantum anharmonic oscillator, whose Hamiltonian is given 
in ([29)1. Dyson's argument would go as follows: If the coupling constant g is continued in 
the complex-g plane to —g, then the potential is no longer bounded below, so the resulting 
theory has no ground state. Thus, the ground-state energy E (g) has an abrupt transition 
at g — 0. If we represent E (g) as a series in powers of g, this series must have a zero radius 
of convergence because E (g) has a singularity at g = in the complex-coupling-constant 
plane. Hence, the perturbation series must diverge for all g ^ 0. While the perturbation 
series does indeed diverge, this heuristic argument is flawed because the spectrum of the 
Hamiltonian (|32|) that is obtained remains ambiguous until the boundary conditions that 
the wave functions must satisfy are specified. The spectrum depends crucially on how this 
Hamiltonian with a negative coupling constant is obtained. 

There are two ways to obtain the Hamiltonian (|32j) . First, one can substitute g = \g\e t0 
into the Hamiltonian (|29|) and rotate from 9 = to 9 = it. Under this rotation, the ground- 
state energy E (g) becomes complex. Evidently, E (g) is real and positive when g > and 
complex when g < 34]. Second, one can obtain ()32|) as a limit of the Hamiltonian 

H = l - V 2 + ~mV + - A gx\ixY {g > 0) (35) 

as e : — > 2. The spectrum of this Hamiltonian is real, positive, and discrete. The spectrum 
of the limiting Hamiltonian ()32|) obtained in this manner is similar in structure to that of 
the Hamiltonian in (}2*9*|) . 

How can the Hamiltonian ()32|) possess two such astonishingly different spectra? The 
answer lies in the boundary conditions satisfied by the wave functions 0„(x). In the first 
case, in which 9 = argg is rotated in the complex-^ plane from to tt, ip n {x) vanishes in the 
complex-x plane as \x\ — > oo inside the wedges — 7r/3 < argx < and — 4n/3 < argx < — n. 
In the second case, in which the exponent e ranges from to 2, <p n {x) vanishes in the 
complex-x plane as \x\ — > oo inside the wedges — 7r/3 < argx < and — tc < argx < — 27r/3. 
In this second case the boundary conditions hold in wedges that are symmetric with respect 
to the imaginary axis; these boundary conditions enforce the VT symmetry of H and are 
responsible for the reality of the energy spectrum. 

Apart from the spectra, there is another striking difference between the two theories cor- 
responding to H in (}3*2*|) . The one-point Green's function G\(g) is defined as the expectation 
value of the operator x in the ground-state wave function (po(x), 

Gi{g) = (0|x|0)/(0|0) = J dxx^l{x) / J dxip%(x), (36) 

where C is a contour that lies in the asymptotic wedges described above. The value of 
Gi(g) for H in (J3~2"j) depends on the limiting process by which we obtain H. If we substitute 
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g = g e l6 into the Hamiltonian ()29|) and rotate from 6 = to 6 = ir, we find by an elementary 
symmetry argument that G\(g) = for all g on the semicircle in the complex-^ plane. Thus, 
this rotation in the complex-^ plane preserves parity symmetry (x — > —x). However, if we 
define H in (|32j) by using the Hamiltonian in (|35|) and by allowing e to range from to 2, 
we find that Gi(g) ^ 0. Indeed, G\(g) ^ for all values of e > 0. Thus, in this theory VT 
symmetry (reflection about the imaginary axis, x — > —x*) is preserved, but parity symmetry 
is permanently broken. We believe that this means that one might be able to describe the 
dynamics of the Higgs sector by using a — gip 4 quantum field theory. 
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